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FOREWORD 


This  report,  L&T  No,  636,  la  the  final  report  under  Contract  1)00421- 
76-C-0052  with  the  Patuxent  River  Naval  Air  Test  Center. 

The  researches  were  carried  out  by  J.D.  Robb  and  Dr,  T.  Chen  with 
major  assistance  from  Wm.  Walker  of  Patuxent  River  Naval  Air  Test  Center 
who  was  also  the  technical  representative  for  the  Patuxent  River  Naval 
Air  Station. 
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ABSTRACT 


Studies  have  been  carried  out  to  detennine  the  hot  spot  temperatures  | 

of  s\;ept  laboratory  lightning  discharges  using  an  infra  red  scanning  camera,  j 

modified  to  provide  higher  time  resolution  with  the  LTRI  St.  Paul  wind  | 

tunnel  used  to  sweep  the  stroke.  The  data  was  compared  with  non  swept  (sta-  ] 

tionary  arc)  data.  Somewhat  higher  energies  were  used  in  the  test  discharges  | 

corresponding  to  the  new  waveforms  (SAJJ-Task  F)  including  a 100,000  restrike  | 

and  a 15  coulomb  intermediate  component  superimposed  on  two  hundred  amperes  5 

of  continuing  current  typical  of  natural  lightning.  The  results  indicate  j 

that  with  the  thinner  new  materials,  a hot  spot  Ignition  hazard  definitely 
exists  for  the  lower  conductivity  materials,  titanium  and  stainless  steel  * 

but  less  than  had  been  previously  estimated.  Yet  to  be  determined  is  the  ' 

complete  envelope  of  hot  spot  time  duration  versus  temperature  required  for 
Ignition. 

i 

i 


i 


1 . 0 INTRODUCTION  AITO  SUMMARY 


I 


One  of  the  areas  of  aa  aircraft  most  Inherently  vulnerable  to  lightning 
strikes  Is  the  fuel  system  and  much  effort  has  been  directed  over  a period  of 
many  years  to  assure  the  safety  of  new  systems  through  severe  test  standards. 
With  the  development  of  advanced  aircraft  having  new  and  thinner  skin  materials, 
the  question  of  possible  heating  of  Integral  fuel  tank  skins  to  the  spontaneous 
Ignition  temperature  of  the  fuel  becomes  of  renewed  Interest. 

A number  of  earlier  experimental  studies  using  only  stationary  teat  dis- 
charges (without  wlndstream)  have  Indicated  that  aluminum  could  probably  not 
cause  hot  spot  Ignition  because  of  Its  excellent  thermal  and  electrical  con- 
ductivity unless  puncture  occurred  and  this  was  confirmed  by  theoretical  studies. 
Of  concern,  however,  are  the  newer  materials  Including  titanium  and  stainless 
steel  with  their  much  lower  thermal  and  electrical  conductivities. 

The  studies  were  carried  out  using  an  Infra  red  scanning  camera,  modi- 
fied to  provide  higher  time  resolution  and  with  the  LTRI  St.  Paul  wind  tunnel 
used  to  sweep  the  stroke.  The  data  was  compared  with  non  swept  (stationary 
arc)  data.  Somewhat  higher  energies  were  used  In  the  test  discharges  corre- 
sponding to  the  new  waveforms  (SA£-Task  F)^  including  a 100,000  restrlke  and 
a 15  coulomb  Intermediate  component  superimposed  on  two  hundred  amperes  of 
continuing  current  typical  of  natural  lightning.  The  results  indicate  that 
with  the  thinner  new  materials,  a hot  spot  Ignition  hazard  definitely  exists 
for  the  lower  conductivity  materials,  titanium  and  stainless  steel  but  less 
than  had  been  previously  estimated.  Yet  to  be  determined  Is  the  complete 
envelope  of  hot  spot  time  duration  versus  temperature  required  for  Ignition. 

2.0  STROKE  HEATING  AID  COOTACT  MECHANISMS 

Lightning  stroke  contact  with  integral  fuel  tank  skins  Is  established 
only  In  swept  strike  zones  where  the  lightning  sweeps  over  the  skin  as  fuel 
tanks  are  not  or  should  not  be  located  In  direct  strike  zones.  The  mech- 
anisrou  of  stroke  contact  heating  are  complex.  Involving  plasma  contact  with 
the  skin  materials.  Ion  bombardment  and  ohmic  heating  of  the  skin  modified 
by  heat  transfer  with  phase  change  and  wlndstream  cooling. 

Earlier  researches  have  determined  the  charge  transfer  magnitudes  as 
a function  of  the  time  required  for  puncture  and  also  required  for  hot  spot 
Ignition  using  stationary  (non  swept)  arcs  (1,2,3).  The  minimum  energy  and 
cliarge  transfers  occurred  at  discharge  times  of  about  10  to  20  milliseconds. 

These  earlier  researches  had  also  shown  typical  stroke  hang-on  times  of  two 
to  five  milliseconds  for  titanium  and  aluminum  respectively  and  even  shorter 
times  for  stainless  steel  (4). 

This  data  along  with  earlier  theoretical  and  experimental  studies  of 
time  temperature  profiles  (5)  have  Indicated  that  without  puncture,  hot 
spot  ignition  with  aluminum  fuel  tank  skins  is  of  very  low  probability. 
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with  the  trend  toward  thinner  skins  and  the  newer  laaterlals,  the  question 
of  hot  spot  Ignition  again  becomes  of  Increased  Interest  and  an  experimental 
program  has  been  carried  out  to  determine  the  time  temperature  waveforms. 


3.0  EXPERIMENTAL  ARRANGEMENT 


The  simulated  lightning  discharge  arcs  were  swept  over  the  test  sample 
using  a wind  tunnel  with  the  test  sample  held  In  a constant  cross  section 
tliroat  area  as  shown  In  Figure  1.  A Eames  Infra  red  camera  was  located  In 
an  aluminum  cabinet  to  provide  shielding  for  the  electronics  with  an  opening 
only  for  the  Infra  red  optical  path.  Four  artificial  lightning  generators 
were  used  to  supply  the  Initial  continuing  current  and  the  combined  high 
current  and  intermediate  current  restrlke.  The  restrlke  was  triggered  by 
a high  current  rate  of  rise  discharge  from  a separate  generator.  The  wave- 
form Is  shown  In  Figure  2 and  a schematic  diagram  of  the  generators  and 
their  connections  Is  shown  In  more  detail  In  Figure  3. 


Approximately  three  test  runs  were  made  for  each  test  sample  at  each 
of  three  levels  of  charge  transfer,  3,  9 and  15  coulombs.  The  maximum  temp- 
erature recorded  was  taken  as  the  measured  value. 


The  test  sample  was  masked  with  tape  to  provide  a narrow  test  strip  of 
bare  metal  along  the  direction  of  air  flow  and  assure  that  the  hot  spots  all 
fell  within  the  narrow  field  of  view  of  the  Infra  red  camera.  To  Improve 
the  time  resolution,  the  Barnes  Infra  red  camera  was  used  with  the  vertical 
scan  disabled  and  with  a narrow  test  sample.  The  IR  scan  was  swept  only 
along  a narrow  strip  of  the  test  sample  at  a rate  of  about  2 milliseconds 
per  scan.  The  data  was  recorded  on  a Hewlett-Packard  Type  3960C  Instrumenta- 
tion tape  recorder  with  a bandpass  at  15  Inches  per  second  of  50  to  60,000 
Hertz.  The  Barnes  camera  had  a resolution  of  10  nllllradlans  equivalent  to 
120  kilohertz  in  the  frequency  domain  but  the  tape  recorder  limited  the  re- 
solution to  about  20  mllllradlans  (60  kilohertz)  corresponding  to  about  5/8 
Inch  of  exposed  width  of  test  panel.  The  theoretical  calculations  Indicated 
a backside  temperature  spot  width  at  peak  temperature  of  slightly  less  than 
5/8th  Inch.  Optical  alignment  of  the  camera  with  the  test  object,  which  was 
critical  to  the  test  result  accuracy,  was  accomplished  by  placing  a heated 
wire  along  the  arc  path  and  observing  the  Barnes  output  with  an  oscilloscope 
until  the  output  was  even  along  the  entire  scope  trace. 


The  Barnes  camera  output  was  calibrated  several  times  during  the  program 
using  a laboratory  standard  "black  body"  temperature  source  with  an  appro- 
priate aperture.  The  output  was  recorded  on  one  channel  of  the  tape  recorder 
and  a synchronization  pulse  was  recorded  on  an  adjacent  channel. 


The  output  had  sufficient  resolution  to  show  clear!' 
rise  at  major  are  hang-on  points  as  well  as  the  spot  wldti 


temperature 
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HIGH  CURRENT  RATE  OF  RISE  COMPONENT 


TIME  - MILLISECONDS 


Figure  2.  Test  Current  Waveform  for  Hot  Spot  Temperature  Measurements. 
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4.0  ANALYSIS  AITO  EXPERIMENTAL  RESULTS 

The  theoretical  aiialysls  of  the  Internal  hot  spot  temperatures  are  com- 
plex because  of  the  different  heating  modes;  joule  heating  and  Ion  bombardment 
and  the  phase  changes  - solid  to  liquid  and  liquid  to  vapor.  In  earlier  LTRI 
theoretical  studies  (Appendix  I and  II) . the  hot  spot  time  temperature  profile 
of  aluminum  was  calculated.  In  these  calculations,  a constant  temperature 
heat  Input  was  assumed  because  of  the  thermal  reservoir  formed  by  a molten 
pool  of  aluminum  and  these  calculations  penalt  a comparison  of  the  new  mate- 
rials coming  Into  use  for  advanced  aircraft.  As  Illustrated  In  Figure  4 from 
L&T  Report  333,  the  simplified  calculation  shows  a faster  rise  and  fall  time 
and  a lower  peak  temperature  than  was  measured  using  a thermocouple  with  the 
stationary  arc  on  aluminum.  Apparently  the  molten  metal  slows  the  heat  rise 
by  holding  the  heat  Input  at  a temperature  just  above  the  melting  temperature 
(659 *C  for  aluminum)  which  Is  well  below  the  arc  temperature  by  at  least  an 
order  of  magnitude.  The  molten  metal  apparently  also  slows  the  temperature 
decay  as  It  solidifies. 

Tests  were  thus  carried  out  using  swept  strokes  on  the  new  skin  mate- 
rials as  outlined  In  Section  2.  In  the  data  reduction  and  analysis  phase 
of  the  present  program,  the  experimental  results  stored  In  the  recording 
tape  were  reviewed  In  two  basic  forms,  one  In  whlcii  the  Individual  traces 
were  shown  as  Illustrated  In  Figure  5 to  permit  examination  of  the  details 
of  hang-on  and  temperature  change  and  one  in  which  only  the  envelope  of  the 
peak  amplitudes  were  displayed  for  determining  the  peak  temperature  and, 
equally  Important,  the  rise  and  fall  times  as  Illustrated  In  Figure  6. 

In  Figure  5 are  shown  two  hang-on  points,  one  solid  and  one  In  which 
the  arc  Is  moving  about  considerably,  i'uch  higher  temperatures  are  recorded 
for  the  point  at  which  the  arc  Is  solidly  attached  as  might  be  expected. 
Generally  the  point  with  the  shortest  total  arc  path  length  nay  be  expected 
to  hang  the  most  solidly. 

In  Figure  6 are  Illustrated  the  great  differences  In  rise  and  fall  times 
and  peak  amplitudes  for  two  different  materials,  aluminum  above  and  stainless 
steel  below.  The  much  faster  rise  and  fall  times  and  lower  amplitudes  of 
the  aluminum  Is  Illustrated.  As  the  minimum  Ignition  temperature  Is  a func- 
tion of  the  time  duration,  or  In  effect  the  total  heat  added,  the  much  lower 
probability  of  fuel  Ignition  with  the  aluminum  Is  Illustrated.  This  con- 
firms the  earlier  experimental  and  theoretical  studies  of  aluminum. 

The  wind  tunnel  tests  have  shown  a great  reduction  in  hot  spot  temper- 
ature with  sweeping.  A comparison  between  the  swept  and  stationary  arc  (for 
a few  metals  run  during  the  present  program  for  direct  comparison)  show  much 
lower  temperatures  by  a factor  of  four  to  six  and  simple  calculations  do  not 
indicate  that  simple  wlndstream  cooling  of  the  skin  would  explain  the  diff- 
erence. Greater  evaporation  of  the  molten  metal  and  arc  movement  might. 

The  measured  values  for  the  few  samples  In  which  both  stationary  and  swept 
tests  were  made  are  presented  In  Table  I. 
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TEMPERATURE 


Figure  A.  Experimental  and  Theoretical  Tine  Temperature  Curves  for 
1/3  Inch  Aluminum  from  Earlier  Studies  (3). 
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HORIZOHTAL  SCALE  - n.l25  SEC/DIV,  VERTICAL  - 314“C/  OIV 


Figure  6.  OscllloRramn  of  Peak  Temperature  Envelopes  Rhov  Great 
Variation  in  Rise  and  Decay  of  0.060  Aluminum  (above.) 
and  Stainless  Steel  (below) . 


TABLE  I 


I 
I 
I 

I Hot  Spot  Metal  Temperatures  With  Comparison  of  Swept  and  Stationary  Arcs 


Aluminum 

Titanium 

Stainless  Steel 

t - 0.060 

t - 0.090 

t - 0.060 

Stationary 

330*C 

408*C 

942*C 

Swept 

59*C 

94 ’C 

157*C 

Ratio 

5.6 

4.3 

6 

A comparison  of  the  temperatures  between  the  LTIU  theoretical  calcu- 
lations and  the  measured  stationary  and  swept  arcs  using  the  infra  red  cam- 
era is  shown  in  Figure  7.  The  figure  illustrates  that  the  swept  stroke  hot 
spot  profiles  lie  well  below  the  LTRI  theoretical  curves  and  much  below  the 
measured  stationary  arc  temperatures. 

A final  summary  of  the  temperatures  as  a function  of  restrlke  charge 
transfer  is  presented  in  Table  II,  As  may  be  seen,  with  swept  strokes,  the 
boron  and  graphite  have  the  liighest  temperatures  followed  by  titanium  and 
stainless  steel  lii  the  middle  range  of  thickness.  It  should  again  be  em- 
phasized that  the  hang-on  times  of  the  materials  substantially  affect  the 
results,  particularly  the  very  low  hang-on  times  of  stainless  steel  which 
greatly  reduces  the  relative  hot  spot  temperature. 

5.0  CONCLUSIONS  AIH)  RECO!CIEiroATIOHS 

Swept  stroke  temperatures  from  hot  spot  heating  by  lightning  current 
arcs  are  substantially  below  the  temperatures  recorded  with  stationary  arcs 
which  have  been  used  in  all  previous  studies  and  on  which  most  fuel  ignition 
hazard  evaluations  have  been  made. 

This  data  permits  a major  reevsluatlon  of  ignition  hazards  with  the 
newer  lower  conductivity  materials  such  as  titanium  and  stainless  steel. 
\fliat  remains  to  be  determli\ed,  however,  is  the  complete  envelope  of  time 
temperature  required  for  fuel  vapor  ignition  and  this  is  recommended  for 
an  additional  study.  This  would  permit  determination  of  safety  margins 
between  the  skin  temperatures  and  the  fuel  vapor  ignition  time  temperature 
profiles. 

Also  needing  to  be  evaluated  are  the  current  time  envelopes  required 
for  skin  puncture  with  swept  strokes  and  an  evaluation  of  hot  spot  and 
puncture  current  time  envelopes  for  painted  skins.  These  are  also  re- 
commended for  an  additional  study  as  most  aircraft  now  require  painting 
for  corrosion  protection. 
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TABLE  II 

Maximum  Recorded  Hot  Spot  Temperaturea 
With  Corresponding  Time  Durations 


est  Sample 


Temi 


C.  Time  Duration  - Seconds 
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1,0  INTRODUCTION  AND  SUMMARY 


Studies  have  been  carried  out  to  determine  theoretical  Integral  fuel 
tank  skin  temperatures  when  subjected  to  lightning  swept  stroke  restrikes. 
The  studies  were  made  for  the  new  fuel  tank  skin  materials  coming  into  use 
for  advanced  aircraft  including  titanium,  stainless  steel  and  graphite  and 
boron  epoxy.  The  new  materials  show  much  higher  peak  temperatures  than 
aluminum  used  heretofore  for  nearly  all  Integral  fuel  tank  walls.  The 
studies,  along  with  current  experimental  measurements  of  fuel  tank  skin 
temperatures,  permit  a better  evaluation  of  hot  spot  ignition  hazards. 


2.0  THEORETICAL  BACKGROUND 

The  studies  are  based  on  earlier  theoretical  work  (Appendix  I)  in  the 
solution  of  the  partial  differential  equation  of  heat  flow  for  a quantity  of 
heat  at  the  melting  temperature  suddenly  added  to  a finite  volume  in  a thin 
plate.  The  basic  assumption  is  that  as  long  as  the  skin  has  not  melted 
through  to  the  inner  surface  that  no  puncture  has  occurred,  the  inner  sur- 
face must  remain  below  the  melting  temperature.  The  solution  is  carried 
out  in  rectangular  coordinates  for  simplicity  but  could  be  converted  to 
cylindrical  coordinates  which  would  not  affect  the  results  substantially. 
This  solution  does  not  recognize  the  complex  mechanisms  Involved  In  the  heat 
transfer  process  through  Joule  heating  and  Ion  bombardment  or  In  the  phase 
changes  from  solid  to  liquid  and  liquid  to  vapor.  These  additional  aspects 
considerably  complicate  the  solutions  and  until  they  are  better  understood 
do  not  contribute  significantly  to  greater  confidence  In  the  results  than 
does  the  simple  classical  heat  flow  solution.  The  simpler  solutions  have 


therefore  been  used  to  permit  ready  comparisons  between  the  new  materials. 
This  data  permits  some  Interpolation  of  new  experimental  data  on  hot  spot 
temperatures  for  both  stationary  and  swept  stroke  arcs, 

3.0  ANALYSIS 

The  results  of  the  calculations  are  presented  In  Figures  1 through  A 
and  are  presented  In  both  log  linear  and  linear  coordinates.  On  each  graph 
Is  also  plotted  an  estimated  minimum  Ignition  tlme~temparature  curve  based 
on  the  few  presently  available  data  points  with  straight  line  Interpolation 
between  the  extremities.  The  plots  are  presented  for  several  materials  of 
the  same  thickness  to  permit  comparison  of  temperature  profiles. 

Of  interest  Is  the  short  duration  of  the  temperature  pulse  of  titanium. 
As  the  total  time  termperature  Integral  determines  the  Ignition  probability, 
titanium  would  seem  to  be  superior  to  stainless  steel  because  of  the  more 
rapid  temperature  drop.  It  should  be  noted  that  In  experimental  studies, 
stainless  steel  has  been  shown  to  have,  however,  a very  short  hang-on  time 
(Appendix  II)  which  balances  Its  longer  duration  hot  spot. 

In  tests  of  0.012"  (12/1000)  nickel  plated  stainless  steel,  no  punc- 
ture occurred  from  100,000  ampere  discharges  with  Intermediate  currents  be- 
tween five  and  ten  coulombs.  This  was  attributed  to  the  low  hang  on  time 
and  arc  spread  characteristics.  There  were  about  100  pit  marks  per  foot 
spread  both  laterally  and  longitudinally  as  shoim  In  Appendix  II.  Incid- 
entally It  should  be  noted  that  the  tubing  was  pressurized  during  the  tests 
to  2650  psl.  Boron  and  graphite,  because  of  their  much  higher  malting 
temperatures,  show  significantly  higher  hot  spot  temperatures  than  titanium 


and  stainless  steel 


From  these  calculations,  an  estimate  of  the  relative  capacities  of  the 


I 
I 

^ various  materials  to  withstand  lightning  discharge  currents  is  shown.  It 

I should  be  noted  that  as  with  the  stainless  steel  swept  stroke  testa  of 

Appendix  II,  a critical  factor  is  sweeping  which  greatly  reduces  the  hot 
I spot  te.iperature.  This  is  also  Indicated  in  Appendix  III,  a paper  on  some 

experimental  measurements  of  inner  skin  temperatures  with  swept  strokes. 

4.0  CONCLUSIONS 

Theoretical  studies  of  fuel  tank  skin  material  hot  spot  temperatures 
caused  by  lightning  restrike  currents  have  shown  the  much  higher  tempera- 
tures reached  by  titanium  and  stainless  steel,  but  also  the  much  shorter 
duration  of  the  titanium  temperature  pulse  as  compared  to  stainless  steel. 
Even  higher  temperatures  are  Indicated  for  boron  and  graphite  epo^. 
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Figure  1.  Calculated  Hot  Spot  Temperature  for  Five  Materials 
of  0.12  Thickness. 
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Figur«  A,  Calculated  Hot  Spot  Temperature  for  Five  Materials 
of  0.03  Inch  Thickness. 
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APPENDIX  II 


HEAT  FLOW  IN  A PLATE  DUE  TO  LIGHTNING  STROKE 


Introduction 

When  a lightning  stroke  hits  a metal  plate,  such  as  a structural 
raomher  of  an  airplane,  it  may  burn  through  the  plate  or  it  may  burn  out 
only  a small  pool  of  metal,  causing  a pitting  of  the  plate.  Of  course, 
the  worst  damage  is  caused  when  a hole  is  burned  through  the  plate.  The 
material  on  the  back  side  is  then  open  to  damage  by  the  arc;  and,  if  a 
member  such  as  a gas  tank  is  involved,  explosion  may  result. 

Even  if  the  plate  is  not  burned  through,  there  may  be  a sufficiently 
great  release  of  heat  energy  in  the  plate  to  lead  to  severe  damage.  The 
question  is  whether  a large  enougli  temperature  rise  may  occur  on  the 
inner  surface  of  the  wall  of  a gas  tank  to  Initiate  an  explosion  in  the 
tank.  Tnis  problem  in  heat  conduction  is  theoretically  analyzed  herein. 


Tlieory 

Consider  a meral  plate  in  the  fonn  of  a square  of  side  a and  of 
tliickness  b,  with  a » b.  Take  a rectangular  set  of  coordinate  axes 
with  origin  at  the  center  point  of  one  side  of  the  plate,  referred  to  as 
the  inside  surface.  Orient  the  x, y-axes  parallel  to  the  sides  of  the 
piate,  and  the  z-axis  normal  to  the  inside  surface.  T}ie  outside  surface 
of  the  plate  then  corresponds  to  the  surface  z = b.  Tlie  lightning  stroke 
will  be  supposed  to  strike  the  outside  surface  of  the  plate  at  the  point 
(0,0,b),  directly  opposite  the  origin  of  coordinates  0 of  the  coordinate 
system,  and  to  release  there  a quimtity  of  heat  Q at  time  t = 0.  The 
problem  is  to  evaluate  the  temperature  distribution  on  the  inner  surface 
of  the  plate  as  a function  of  time;  in  particular,  the  temperature  at 
the  point  0,  which  will  be  the  hottest  point  on  the  inner  surface. 

Since  this  is  a damage  problem,  it  is  best  to  set  the  conditions  on 
the  analysis  so  that  the  temperatxire  rise  on  the  inner  surface  will  be 
overestimated  rather  than  underestimated.  Therefore,  the  analysis  is 
simplified  by  supposing  that  no  heat  is  lost  from  the  plate  either  by 
radiation  or  by  conduction  to  the  air. 

The  heat-conduction  equation 

(.2»2  - ^)  T - 0 (1) 


I 


is  to  be  solved  in  the  region 


-a  ^ a 

2 - - 2 


-a  a 

— < y < — 

2 - - 2 


0 < z < b 


with  the  boundary  conditions 


= 0 

at 

II 

X 

= 0 

at 

y = i- 
2 

= 0 

ab 

z = 0,b 

The  following  functions  are  special  solutions  of  the  heat-conduction 
equation  (l)  that  satisfy  the  boundary  conditions  (3): 

+ f)  [©(y  ^ f)  exp(-Xi^t)  (4) 

where  2,  m,  and  n are  non-negative  integers  and 


^2mn  - 


fe'  ^ (?/  * (.f/] 


Next,  these  particular  solutions  are  combined  in  such  a manner  that 
the  initial  conditions  representing  the  heat  source  produced  on  the  out- 
side of  the  plate  by  the  lightning  stroke  are  satisfied.  For  this  pur- 
pose the  general  solution  for  the  temperature  at  any  point  in  the  plate 
is  formed: 


T(x,y,z,t)  = 


Azmnf 


Z n^O 

Let  the  temperature  distribution  at  t = 0 be  known,  so  that 

T(x,y,z,0)  = u(x,y,z) 
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is  a known  function  of  the  coordinates  x,  y, z in  the  plate.  Combining 
equations  (6)  and  (v)  gives 


u(x,y,z)  = 


^Imn 


In  order  to  express  the  coefficients  in  this  Fourier  series  in  terms 
of  the  function  u(x,y, z),  the  symbol  M[F(x, y, z)]  is  defined  as  the 
mean  value  of  the  function  F(x, y,  z)  over  the  volume  of  the  plate.  That 


M[F(x,y,z)]  = 


/•a/2  ra/2  rb 

J.a/2  J-a/z  Jo 


F(x, y, z)  dz  dy  dx 


^ J-a/Z  J-a/z 

where  a^b  is  the  volume  of  the  plate. 


Using  equation  (S)  and  the  usual  methods  of  calculating  the  coeffi- 
cients in  a Fourier  series  gives 


^ ^ M{u(x,y,z)  cos  ['[M)(x  . |)]  cos  | (M)[y  . |)|  cos  | (M)z|} 

M {cos"  [(^)(x  » I)]  00.2  [[f){y  . I)]  00.2 

(10) 

The  evaluation  of  the  initial  temperature  distribution  u(x, y, z) 
requires  simulating  the  heating  effect  of  the  lightning  stroke  in  some 
suitable  fashion.  Suppose  that  at  t = 0 an  amount  of  heat  Q is  lib- 
erated in  a small  block  of  the  plate  defined  by  the  conditions 


< x<  ^ 
2 2 - — 2 2 


-=<— <y<-<-  / ui; 

2 2 2 2 

b-P<.z<b  ^ 

T/ie  heat  is  thus  liberated  in  the  small  block  of  dimensions  a, a, 3 sit- 
uated on  the  outer  surface  of  the  plate,  just  opposite  the  point  0 on  the 
plate.  The  material  in  this  block  will  be  raised  to  the  uniform 
temperature 
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0 = -V 

cpa^3 

where  c is  the  specific  heat  and  p is  the  density  of  the  metal.  The 
function  u(x,y,z)  will  he  eclual  to  0 in  this  block  and  equal  to  zero 
outside  it. 

The  coefficients  from  equation  (lO)  can  now  be  evaluated  as 

follows: 


li. 

•a/2 

COB 

-a.!'?. 

•f)hl 

\i) 

ra/2 

COS 

-'V2  . 

It) 

(-1) 

dy  1 

r 

i 1 COS 

b Jb-ti 

It)- 

dz 

J t - ~ 

i(l  + 6, 
2' 

_l(l  - ^ 

,o\ 

^n,o) 

(13) 


For  convenience  of  notation  the  Kronecker  symbols  are  used: 

fl  if  1 = 0 
6 = < 

\o  If  2 i(  0 

Carrying  out  the  integrations  indicated  in  equation  (13)  yields 

t [(?)(^  ^ t)]  ^ H'*  fe)(x)]  ■ 

J-CLlZ  '■ 


(-) 


sin 

il^) 

(t) 

if 

1 is  odd 

,1/2 

(r) 

- if 

Ina 

a 

2a 

is  even 


(15) 
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■ (f^)  =“["  (fir)_ 

“ <->t;l)  (r)  {^) 

^ \b/  nitP 


This  yields  the  general  formula 


^Zmn  ~ ® 


[|(i-h,o)][l<-Vo)][|(i-v„o)]  ir)  (#)  (S^) 


Here  Z and  m have  only  even  values,  since  by  equation  (lb)  the  coef- 
ficient will  vanish  if  I or  m is  odd. 

The  final  solution  of  equation  (6)  for  the  temperature  distribution 
in  the  plate  for  t > 0 can  be  written  as 

T(x,y,z,t)  = — %-  <P(a,a;x,t)<P  (a,a;y,t)  \^(b,3;z,t)  (18) 

cpa^^P 


<p(a,a;x,t)  = 


2—0 7 2 ^ 4^  • • • 


* ^,o)  If 


i-m'] 


i(b,3;z,t) 


(-)"£c<,srWzl  sin  (»«!!) 


n=:0,l,2. 


|h  * 


(20) 
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If  I is  even. 


cos 


(21) 


so  that  the  new  summation  index,  p = Z/2,  can  be  introduced,  which  takes 
all  non-negative  integral  values,  and  formula  (19)  can  be  replaced  by 


<p(a,a;x,t) 


p=0, 1,2,  • • 


+ 8z,o) 


(22) 


The  foraiulas  given  represent  the  exact  solution  of  the  proposed 
problem.  However,  in  view  of  the  nature  of  the  applications  to  be  made 
of  them,  it  is  feasible  to  allow  some  simplifications. 


Hie  size  of  the  plate  is  not  very  material  to  the  flow  of  the  heat 
to  the  Inside  of  the  plate  so  long  as  the  thickness  is  small  compared 
with  the  breadth  of  the  plate.  It  is  therefore  convenient  to  take  the 
limit  a-»«».  In  this  case, 

<p(a;x,t)  = lim  <p(a,a;x,t)  (23) 

a “ 

The  series  in  equation  (22)  goes  over  into  an  integral,  and  on  working 
out  the  process  the  following  solution  is  found: 


<p(a;x,t) 


(24) 


where  ^ = pjta/a. 

This  function  can  be  reduced  to  a more  convenient  form  for  numerical 
computations.  Introduce  the  following  notation  (which  is  veilid  for 
t > 0); 


(25) 


i 


Substitution  into  equation  (24)  yields 


exp  (-u^) 


du 


Wow  define  the  function 


g(n) 


^ 2 f sin_,(^),  (.^2)  ^ 

Jn  ^ 


(26) 


(27) 


Making  use  of  this  abbreviation  gives  from  equation  (26), 


<p(a;x,t) 


(28) 


From  equation  (27), 

AqD 

= i I COS  (2Ti^t)  exp  (-p^)  d+i 

« Jo 

= exp  (-n^)  (29) 

on  making  use  of  formula  508  of  Peirce's  table  of  integrals  (ref.  2). 

It  is  obvious  from  equation  (27)  also  that  g(0)  = 0.  Equation  (29)  is 
integrated  directly  as  a differential  equation: 


s(h)  = 


exp 


(-^2)  dp 


(30) 


Tn is  shows  that  the  function  g(q)  defined  in  equation  (27)  is  just  the 
ordinary  probability  integral  and  thus  can  be  found  tabulated  in  numeri- 
cal form. 


;i 

|i 

!i 
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The  following  formula  for  the  temperature  distribution  from  the  heat 
source  (strictly  speaking,  in  a plate  of  infinite  breadth)  is  finally 
obtained: 

T(x,y,z,t)  = <P(a;x,t)<p(a;y,t)iir(b,3;z,t)  (3l) 

cpa'^P 

with  the  function  (p(a;x,t)  defined  by  equation  (28)  and  the  function 
\l((b,  3;z,t)  defined  by  equation  (20). 


Discussion  of  Formula  (31) 

Formula  (3l)  gives  the  temperature  distribution  throughout  the  plate 
resulting  from  the  heat  source  introduced  by  the  lightning  stroke.  In 
applying  the  result,  only  the  temperature  distribution  over  the  inner  sur- 
face needs  to  be  known.  This  is  obtained  by  setting  z = 0 in  equation 
(3l).  For  simplicity  of  notation,  this  temperature  function  is  written 
as  follows: 


To(x,y;t)  = T(x,y,0;t)  (32) 

The  following  formula  results: 


To(x,y;t) 


cpa^p 


<p(a;  X,  t)  <p  (a;  y,  t)  (b,  3;  0,  t ) 


(33) 


where,  from  equation  (20), 


4r(b,3;0,t)  = ^ 


ri—  0 ^ X ^ • • • 


■'  (.)»  sin  (2^) 


exp 


(34) 


Since  the  temperature  distribution  TQ(x,y;t)  is  expressed  in  equa- 
tion (33)  as  the  product  of  two  types  of  functions,  it  is  convenient  to 
examine  the  nature  of  each  of  these  functions  separately. 


At  the  initial  Instant  t = 0,  conditions  require  that  the  whole  of 
the  inner  surface  of  the  plate  be  at  the  uniform  temperature  taken  to  be 
T = 0.  Therefore,  the  following  must  be  obtained  from  equation  (34) : 


Mr(b,3;0,0)  = ^ 

D 


0 


(35) 
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ITiiis  is  an  alternating  series  of  a type  somewhat  difficult  to  handle, 

since  the  terns  do  not  diminish  very  rapidly  in  magnitude.  There  will  be 
no  attempt  to  prove  rigorously  that  the  sum  of  the  series  is  actually 
zero  as  is  indicated  in  equation  (35). 

I 

It  is  apparent  by  inspection  of  equation  (34)  that  as  t -♦  » the 
series  converges  rapidly  to  the  value 

V(b,3;0,-)  = (36) 

Owing  to  the  exponential  nature  of  the  summands  in  the  series  in  their 
dependence  on  the  time  variable,  the  dominating  term  will  be  the  one 
having  the  smallest  exponent;  that  is,  the  second  member  of  the  series, 
since  the  first  one  (n  = O)  does  not  depend  on  the  time  at  all.  There- 
fore, the  function  given  by  equation  (34)  will  start  from  zero  at  t = 0 
and  will  rise  quickly  to  the  final  value  P/b  practically  like  an  expo- 
nential function  with  the  time  constant 

Expressed  in  physical  terms,  this  function  determines  the  flow  of  heat 
from  the  initial  heat  source  directly  through  the  thickness  of  the  plate. 
The  time  constant  (37)  can  therefore  be  expected  to  be  quite  small  for 
plates  of  ordinary  thickness  such  as  are  used  in  the  construction  of 
aircraft. 

For  an  aluminum  plate, 

Thermal  conductivity  = k = 0.504  cal/ (cm) (sec) (°C) 

Specific  heat  = c = 0.217  cal/(g)(°C) 

Density  = p = 2.70  g/cm^ 

From  this  information,  the  thermal  diffusivity  is 

= — = 0.086  cm^/sec 
cp 

X = /\/ 0. 86  = 0.  93  cm/  //sec 

Taking  b = inch  = 0.318  centimeter,  the  time  constant  is 


0.0119  sec 


(38) 
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The  rise  time  of  the  temperature  on  the  inner  surface  of  the  plate,  di- 
rectly opposite  the  heat  source,  should  thus  be  of  the  order  of  12  milli- 
seconds. This  result  will  not  be  particularly  sensitive  to  the  size  of 
the  initial  heat  source;  that  is,  it  does  not  depend  greatly  on  the  value 
of  p,  since  the  expression  (37)  for  t'  does  not  involve  this  parameter. 

Tlie  temperature  at  the  point  0 should  have  a maximum  value,  which 
is  about 


T 

0,max 


p 

cpa  b 


(39) 


The  initial  rise  of  the  temperature  at  the  point  0 is  determined  by 
the  flow  of  heat  through  the  thickness  of  the  plate,  but  its  ultimate 
decline  is  governed  by  the  transverse  flow  of  heat  along  the  plate.  This 
is  expressed  by  the  functions  (p(a;x,t)  and  (p(a;y,t),  which,  of  course, 
have  the  same  functional  form. 


From  the  initial  conditions,  the  following  is  expected  at 


<P  (a';x, 


if 

if 


|x|  > a/2 
I x I < a/2 


t = 0: 


(40) 


It  is  easy  to  show  from  equation  (28)  that  this  is  the  case,  if  it  is 
noted  from  equation  (30)  that 


gC-n)  = -s(ti) 

£(”)  = 1 

The  fimction  (p(a;x,t)  is  roughly  exponential  in  form  and  decays 
COTT.paratively  slowly  in  relation  to  the  initial  rapid  rise  of  temperature 
at  the  point  0.  It  is  best  shown  in  the  form  of  graphs  drawn  for  specieil 
cases. 

The  analysis  has  assumed  that  the  initial  heat  source  has  a square 
cross  section  and  a depth  3.  In  practice  one  will  have  little  or  no 
control  over  the  exact  shape  of  the  region  in  which  a lightning  stroke 
develops  heat  in  the  plate;  and  on  the  whole  one  will  probably  find  a 
circular  or  roughly  elliptical  spot.  For  points  near  the  center  of  the 
spot  and  points  away  from  the  spot  by  distances  large  compared  with  the 
radius  of  the  spot,  the  exact  shape  of  the  spot  will  be  immaterial.  For 
an  exactly  circular  spot  the  analysis  can  be  made  in  polar  coordinates. 
Thiis  analysis  is  discussed  in  a later  section  for  completeness,  but  its 
use  would  require  numerical  work  with  Bessel  functions,  which  was  not 
considered  justifiable  in  view  of  the  uncertainties  in  the  data. 
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Graphical  Example  | 

In  order  to  show  the  nature  of  the  heat  flow  in  the  plate,  a partic- 
ular case  is  presented  in  graphical  form.  Consider  an  aluminum  plate  with 
the  following  dimensions: 

Thickness  = b = 0.318  cm  (l/8") 

!a  = 0.  5 cm  i 

i 

P = 0.159  cm  (l/l6")  I 

The  source  is  then  assumed  to  be  0. 5 centimeter  square  and  extends  half- 
way through  the  plate.  From  the  data  given  in  the  preceding  section  for 

aluminum,  the  temperature  of  the  source,  for  a given  heat  input  Q,  is  ♦ ' 

initially 

0 = 10.3  Q (42) 

where  Q is  expressed  in  joules,  and  9 is  mne  temperature  rise 
above  room  temperature  in  degrees  centigrade. 

Temper at\ire  at  point  0.  - The  temperature  at  the  point  0,  which  is 
on  the  inside  of  the  plate  just  opposite  the  center  of  the  source,  is 
considered  first.  Making  use  of  equation  (33)  gives 

Tq  = To(0,0;t) 

= e|<P(a;0,t)|  2x,i,(b,|;0,t)  (43) 

The  function  (b,b/2;0,t) , which  determines  the  flow  of  heat  through  the 
thickness  of  the  plate,  is  plotted  in  figure  26.  It  starts  from  zero  and 
rises  to  0.5,  since  the  flow  of  heat  directly  through  the  plate  would 
double  the  amount  of  heated  metal  and  so  would  lower  the  temperature 
by  a factor  0.5.  This  function  is  plotted  on  a universal  time  scale 
as  a function  of  t/t',  where  t'  is  defined  by  equation  (37). 

The  function  <P(0.5;0,t),  which  determines  the  flow  of  heat  away  from 
the  point  0 along  the  plate,  is  plotted  in  figure  27.  This  graph  starts 
at  imity  at  t = 0 and  diminishes  comparatively  slowly  to  zero. 

The  composite  result,  giving  the  temperature  at  the  point  0,  is 
plotted  as  curve  A in  figure  28,  which  shows  that  the  temperature  at  0 
rises  to  about  6/4  as  its  triaximum  value  in  about  20  milliseconds  and 
then  falls  rather  steeply;  in  l/lO  second  it  is  down  to  6/l0.  This 
point  will  obviously  be  the  hottest  on  the  inside  of  the  plate,  so  that 
the  rapidity  with  which  it  cools  off  will  be  an  important  criterion 
governing  the  firing  of  an  explosive  gas  mixture  that  contacts  the  sur- 
face here. 
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Temperat\ire  at  a neighboring  point.  - As  an  indication  of  the  tem- 
peratures reached  on  the  inside  of  the  plate  near  the  source  point,  the 
temperature-time  curve  has  heen  plotted  for  a point  0.5  centimeter  from 
the  point  0.  The  curve  is  given  as  curve  D of  figure  28.  The  tempera- 
ture at  this  point  rises  slowly  to  only  about  25  percent  of  tiie  maximum 
temperature  at  0 and  then  falls  slowly.  The  flow  of  heat  along  the 
aluminum  plate  is  so  rapid  that  only  the  points  quite  near  the  initial 
source  are  heated  to  any  great  extent. 


Effect  of  Continuous  Source 

It  has  been  assumed  in  the  previous  calculations  that  the  lieat  source 
is  established  instantaneously  at  time  t = 0 and  that  only  the  tempera- 
ture distribution  from  this  origin  is  significant.  In  practice,  the 
application  of  the  heat  will  not  be  so  instantaneous,  but  the  source  may 
be  applied  for  some  time  and  may  vary  from  instant  to  instant  in  magni- 
tude. Once  the  problem  of  finding  the  temperature  distribution  from  an 
instantaneous  source  has  been  solved,  it  is  possible  to  write  the  method 
oi"  finding  it  from  a variable  source,  taking  advantage  of  the  linearity 
of  the  differential  equation  of  heat  conduction  and  of  the  boundary 
conditions. 

First  the  notation  in  which  the  result  has  been  expressed  will  be 
revised.  Instead  of  using  the  particular  instant  t = 0 as  the  time  of 
application  of  the  source,  this  instant  is  indicated  as  t^.  Also, 
suppose  that  the  heat  is  supplied  in  an  infinitesimal  time  interval  dtj^, 
such  that  Q = q(t^)dt2^.  Then,  from  equation  (3l)  the  temperature  dis- 
tribution following  from  this  source  at  tines  t > t-j^  would  be  given  by 
the  formxila 

) 

T(x,y,z,t)  = g-  ip(a;x,t-t^)<p(a;y,t-t^)it({b,3jz,t-t^)dt^  (44) 

epa  3 

Clearly,  to  find  the  temperature  distribution  from  a set  of  sources 
operating  in  the  past  it  is  necessary  only  to  sum  (integrate)  expression 
(44)  over  all  the  sources  that  have  been  present.  This  yields  the  final 
formula: 


T(x,y,z,t) 


1 

cpa^3 


q(t^)tp(a;x,t-t^)<p(a;y,t-t^)\y(b,3;z,t-t^)dt^ 


(45) 


It  is  not  practicable  to  evaluate  this  expression,  by  actual  inte- 
gration if  the  source  function  q(tp)  is  very  complex.  The  best  method 
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of  handling  it  is  prohably  by  nvmierical  and  graphical  means,  approximating 
the  source  function  by  a set  of  discrete  sources. 


Treatment  of  Infinite  Plate  in  Cylindrical  Coordinates 

The  mathematical  analysis  of  the  heat  flow  in  a plate  has  been  car- 
ried out  entirely  in  terms  of  Cartesian  coordinates  in  the  earlier  sec- 
tions of  this  report.  This  has  led  to  the  use  of  a heat  source  in  the 
form  of  a small  rectangular  parallelepiped.  The  reader  may  consider  that 
it  would  be  more  sensible  to  use  a heat  source  in  the  form  of  a small 
cylinder.  This  is  certainly  correct  in  principle,  and  the  procedure  used 
has  been  one  of  convenience  only.  In  this  section  the  analysis  is  carried 
through  for  cylindrical  polar  coordinates  with  a cylindrical  shape  for  the 
source  and  solved  in  terms  of  Bessel  functions. 

Using  the  usual  cylindrical  polar  coordinates,  with  origin  at  the 
point  0,  the  heat-conduction  equation  takes  the  form 


(46) 


Only  cylindrically  symmetric  solutions  need  be  considered,  so  that  the 
temperature  depends  only  on  the  distance  from  the  center  of  the  plate 
and  not  on  the  angular  position  around  the  source. 

First,  particular  solutions  of  the  differential  equation  (46)  are 
sought  which  obey  the  boundary  conditions  of  the  problem.  Here  an  in- 
finitely large  plate  is  taken  at  the  start,  so  that  the  boundary  condi- 
tions reduce  to  the  requirement  that  the  solution  be  finite  everywhere, 
and  be  single -valued.  There  is  to  be  no  flow  of  heat  from  the  s'irfaces 
of  the  plate,  so  that  3t/Sz  =0  at  z = 0,b. 

Particular  solutions  satisfying  these  conditions  are  of  the  form 


F(r)  cos  exp  -{^)  tjexp,(-x2t)  (47) 


where  X is  an  arbitrary  real  positive  constant,  and  F(r)  is  a solu- 
tion of  the  differential  equation 


1 

r dr 


0 


(48) 


The  only  solution  of  equation  (48)  that  remains  finite  at  r = 0 (for 
X ^ O)  is 


I 
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F(r)  = constant  X jQ(Xr) 


(49) 


where  Jq  is  the  Be'^sel  function  of  first  kind  of  order  zero. 


To  satisfy  the  initial  conditions  these  particular  solutions  are 
comhined  linearly.  The  members  in  the  vairiable  z must  be  summed  over 
the  non-negative  integer  n = 0,1,2  . . .,  while  the  radial  solutions 
must  be  integrated  over  the  parameter  X.  If  the  initial  heat  source  is 
a small  cylinder  of  radius  a and  depth  P,  into  vhich  an  amount  of  heat 
Q is  deposited  at  t = 0,  the  result  is 


T(r,z,t)  = — (r(a;r,t)\lf(b,P;z,t) 
cpna  P 


(50) 


with 


„ , cos  z sin 

2(^  + ^n,o)  “ 


exp 


(51) 


which  is  identical  with  the  function  defined  in  equation  (20).  The  func- 
tion <p(a;r,t)  is  of  the  form 


<p(a;r,t)  = f g(X)  JQ(Xr)  exp(-x2t)X  dX 

Jo 


(52) 


where  g(X)  irust  be  determined  from  the  initial  i,onditions.  Here  the 
following  is  required; 


<p(a;r,0) 


0 if  r > a 

1 if  r < a 


If 


(p(a;r,0)  = u(r) 


Then,  from  equation  (52), 


(r)  = f e(X)  JQ(Xr)X  dX 
•'O 


(53) 


(54) 


(55) 


The  inversion  of  this  integral  equation  for  g(X)  when  the  left  side  is 
a Known  function  gives 
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s(x) 


u(r)  jQ(xr)r  dr 


(56) 


Making  use  of  the  initial  conditions  (53)  which  define  the  function 
u(r),  in  the  present  problem 


g(x)  = f Jo(Xr)r  dr 

Jo 


This  leads  to  the  formula 


(57) 


<p(a;r,t) 


Jf,(XR)R  dR 


J„(Xr)  exp(-X^t)X  dX 


Jq(XR)  jQ(Xr)  exp(-x2t)X  dX 


R dR 


(58) 


It  would  be  possible  to  make  use  of  these  formulas  for  the  calcula- 
tion of  the  temperature  distribution  in  the  plate,  but  the  work  would  be 
greater  than  by  the  earlier  method,  without  significant  increase  in 
accuracy  of  the  result. 
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